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ABSTRACT 


Quantum phase estimation is a fundamental algorithm in quantum computing for estimating the phase of an 
eigenvalue for a unitary operator. While the output state of quantum phase estimation circuit exactly matches 
to the phase represented as rational number with denominator as powers of 2, the measurement could be read 
in many different states for phases outside this form. In this paper, we address this gap by formulating and 
precisely computing the quantum state of QPE when the phase is not exactly represented as a rational number 
with a denominator in a power of 2. Leveraging mathematical techniques from number theory and 
trigonometry, we derive explicit expressions for the quantum state of QPE in these cases. To validate our 
results, we compare our theoretical predictions with simulations using the Qiskit framework, confirming the 
accuracy of our formulations. Our findings provide insights into the behavior of QPE for phases, expanding 
the understanding and applicability of this fundamental quantum algorithm. 


Keywords: Quantum Phase Estimation, Quantum Fourier Transform, Shor’s algorithm, Probability 


Distribution, Quantum State 
1. INTRODUCTION 


Quantum phase estimation (QPE) stands as a 
foundational algorithm in quantum computing, 
essential for estimating the phase of unitary 
operators with high precision. Originally introduced 
by Kitaev [1], QPE has garnered significant attention 
for its role in various quantum algorithms, including 
Shor's algorithm for factoring large integers 
[2][3][4] and broad range of quantum simulation 
techniques [5][6][7]. Conventionally, QPE is well- 
understood when the phase is represented as a 
rational number with the denominator in a power of 
2, allowing for efficient quantum circuit 
implementations. However, the behavior of QPE for 
phases outside this conventional form remains less 
explored, posing challenges for understanding and 
leveraging the algorithm in practical quantum 
computations. 


In this paper, we address this gap by formulating 
and precisely computing the quantum state of QPE 
when the phase is not exactly represented as a 
rational number with a denominator in a power of 2. 
Representation of output state of QPE in such cases 
is not explicitly formulated in any previous literature. 
Building upon mathematical techniques from 
number theory and trigonometry, our result shows 
that the probability of a state closest to the 
approximation is highest in these cases. Our 


formulations provide a comprehensive 
understanding of QPE's behavior for phases beyond 
the conventional rational representations, offering 
insights into its theoretical foundations and practical 
implications. 


To validate our theoretical results, we employ 
simulations using the Qiskit framework [10], a 
leading platform for quantum computing research 
and development. By comparing our theoretical 
predictions with Qiskit simulations, we confirm the 
accuracy and efficacy of our formulations, 
establishing their utility in practical quantum 
computations. Our findings not only extend the 
understanding of QPE but also pave the way for 
exploring its applications in scenarios where phases 
are not strictly rational with denominators in powers 
of 2. 


2. QUANTUM PHASE ESTIMATION 


Given an n-qubit quantum gate U, Figure 1 
shows a circuit diagram implementing Quantum 
Phase Estimation for gate U. Throughout this paper 
in the circuit diagram, we will use a convention of 
numbering qubit subscripts from bottom to top, 
which means that bottom-most qubit represents the 
highest bit when the measured state is interpreted as 
a binary number. We will name qubits with prefix q 
as in qx; in order to represent the position of the 


aE ERE emmemenimmemimemmmnteseseciee =e ee ee 
3699 


Journal of Theoretical and Applied Information Technology 
30" April 2024. Vol.102. No 8 


SZ 


© Little Lion Scientific 


ISSN: 1992-8645 


qubit in the quantum circuit, where x; represent the 
state of the qubit. In the QPE circuit, the bottom t 
qubits will be measured, and the phase is estimated 
based on the measurement. When Dirac notation 
[11] is used to describe quantum operation, the 
bottom qubit is located to the left of the top qubit. 


HT represent a Hadamard gate, whose 


quantum operation is described with a matrix 
1 1 


4 2 . For an n-qubit quantum gate U whose 
v2 v2 


operation is described with a 2” x 2” matrix U, U* 
represents an n-qubit quantum gate whose operation 
is described with a matrix U*. U gate connected to a 
qubit using a vertical line and a black dot represents 
a controlled-U gate, cUgx,gy, where the control 
qubit is qx; qubit and the target qubits are gys. The 
operation of this (n+/)-qubit gate CUg,,qy is as 
follows. 


for |x;) = |0) 
for |x;) = |1) 


“enabDD)= Eup 


As derived in [12][13] the inverse Quantum 
Fourier Transform circuit, denoted as QFTt , 
transforms the input state to the output state as in the 
following. 
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If the state of gx qubits is measured along the z-axis 
= t = x . 
AS X = Xp Hp_g 1 HyXp E {O11}, A = a: 


phase of eigenvalue of U, which satisfies U]w) = 
e278 ls) for some eigenvector |iy) of U. 


For example, if A= =, the probability 
distribution of measurement of QPE circuit with t=3 
qubits is as in Figure 2. In the case when the 
denominator of A is 2°, the output state is |1) with 
probability = 1, theoretically. 


tA= ~ where the denominator of A is not 
represented as 2°, the probability distribution of 


measurement is as shown in Figure 3. We formulate 
the probability distribution in the following section. 


3. ANALYSIS OF QPE MEASUREMENT 


By the definition of a controlled-U gate, if 
we apply a controlled-U gate, cUgy,qy, to an 
arbitrary quantum state |x;) = a|0) + B|1) and 
ly) = |W), an eigenstate of U, the quantum state 


evolves to the following equation, for an eigenvalue, 
2TtAi 
aia 


CUgx,ay|%i) |b) = al0)lp) + B]L)UlW) 
= «|0)lp) + Bl1)e2"!|y) 
= (al0) + Be?™*|1)) |p) 


Ml > Xo 


|0)S¢- 


Ol - xX) 


—> Xe-1 


ld) Ix) = 12%A) 


Figure I Quantum Phase Estimation circuit used to estimate the phase of an n-qubit quantum gate U 
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The result of this operation changes the 
quantum state of the control qubit when the initial 
state of the target qubit was an eigenstate of U. This 
operation 1s referred as a phase-kickback and used as 
a common technique for measuring operators in 
quantum error correcting codes, such as the surface 
code [14]. The circuit in Figure 1 makes use of the 
phase-kickback. 


Following the equation, 


k k_ 
U* |y) = U* “*Uly) 
ye 1 p7mt iy) 
y2*-2Ye2mhily) 
= y2*-2 ermee ty) 


Y%e2n2*rily) 
ky; 
= e202 Alay) 


|W), an eigenstate of U, is also an eigenstate 


k , ky; 
of U? for eigenvalue e277“, 


The 
Ue a gate is described as in the following. 


operation of the — controlled- 


k 
cue aay eI) 
7 Ixi)ly), 
— k kj; 

|x;)U? ly) = e?™ "| x;) Ly), 


1.000 phase = 1/8 
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Figure 2 Probability Distribution of Measurement in OPE 
Circuit for a Gate with Phase A =< USING 3 


Measurement Qubits. 


for |x;) = |0) 


This implies that he role of the controlled- 


k. ‘ : : 
U*’ is to change the state of control qubit, which is 
one of Xo, QX4,-»,QX--1 qubits in the circuit in 
Figure 1. 


The initial state of the circuit is prepared as 
|0)®*lw) where |W) is an eigenstate of U. 
Hadamard gates H changes the quantum states of 
qx; s |0)>|+)= = (|0)+ |1)) .Thus, the 
quantum state |d) which includes only bottom t¢ 
qubits, @QX;z~-1,+-,9Xo, 1S represented as_ the 
following equation. 


|0) + e2mhi2'*)4) |0) + e2mhi2*)4) 


lp) = —————_ ® ——— ® 
V2 v2 
0) + e2mhi2" 11) 0) Ea e2mhi2° 11) 
; V2 V2 
= 3 > e 2miA(2°~t xp_4+--+2°x9) co 2X5) 
xE{0,1}€ 


using the representation 
1K 1K 
10) a e2mai2 |1) — e2TAi2 |x), 
xeE{0,1} 


The inverse Quantum Fourier Transform 
circuit, denoted as QFT™, transforms the input state 
to the output state and the state is rewritten as in the 


for |x:)=|1) following equation. 


phase = 1/3 
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Figure 3 Probability Distribution of Measurement in OPE 
Circuit for a Gate with Phase A == USING 3 


Measurement Qubits. 
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Xx, (XxX x 
10) +e 211) _ joy +e 2™lat*22) 1) 
7 V2 V2 t-1 t-1 
—2mi (~E5147t524...470) 1 2i0,x 1 2i0 
0) +e Uae a 28/1) co ==] | Dd. zee =5] [+ e2*) 
2 k=0 xp€{0,1} k=0 
1 : —2m1Yt-1X0 : ; ; 
= ra > e 2 Ivt-1) ®@ using the following relation 
Yt-1=0 
il ~2mtiye—p(2%1 +%9) 1+ e78 =1+4cos20+ isin 20 
ion ». © ‘ IYr-2) @ ». = 1+(2cos*6—1) + 2isin®@ cos 6 
Yt-2=0 = 2cos@(cos 6+ isin®) 
1 —2miyo(2°txp¢_44+--+2x14X9) = 2 COS 6 et 
vos > e 2 Ivo) 
= ee eee The probability amplitude of |y),c(y) is 
a > ee... 2 SOS:~C~CSO «Sf Presented as a product of cosine values and 
Ze conyt exponential functions of pure imaginary numbers as 
os in the following. 
Then QFT'|) is represented with c(y), a if 
ey) = =| |G + e7/8k) = | | cos 6; ek! 
QFT"|¢) 2° 
1 k=0 k=0 
ae eo 2mia(2 ~txp_4 4+ +2°X9) QFT* |x) 
xefo1jt Because the square of L2-norm of 
_ probability amplitude indicates the probability of the 
~ cy) ly) state is measured, and L2 norm of exponential 
ye{o,1}¢ 


c(y) as defined in Equation 1 below. c(y) 
is the probability amplitude of measuring |y). 


We define 0; as 


—— 

O28 = Vpgged oe i 2 et G2") 
{2x 

= (2) ane ee 


1 
c(y) ~ 5¢ 


xe{0,1}6 
t-1 


1 
“all 


k=0 xp €{0,1} 


e 2tOkXk 


function of a pure imaginary number is | by Euler's 
formula, probability of measuring |) = |y) is 
analytically computed as a square product of t cosine 
values as in the following equation. 


t-1 


P(ly)) = lei? = | | cos? 4, 


k=0 


For example, in QPE circuit of U gate with 
t=3 qubits for measurement, 


c(y) = cos, e®2' cos 0, e®! cos 0, e%! 


2mi(A(2°-2xp_4++-+20x9)—(ye-12 x9 +--+ 79 (2° xp a +-+X0))/2°) 


(1) 
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where 


0, = (227A — yp 27/23) tt 


1 
(22 > Vane fe TT 


m=0 
= (2'A-— (9 2* +9, 27) / 2°) 0 


(2% - S92 /2 1 


m=0 


= (2° — (Yo a T+ V4 2* + y2 2°) J 2°) in 


0, = 


A 


for U gate with phase A = - , where A # _ 


for any integer k with t=3 qubits for measurement, 
Table 1 shows the probabilities, ||c(y)|| * along with 
65,94, Oe . 


This result matches the graph in Figure 3. 
Because - is close to =, the probability of measuring 


ly) = |011) is highest and the probability decreases 
as the difference between y and 2°A increases. 


Table 1 Probability of Measuring \y) in OPE Circuit for 
a Gate with Phase A = - Using 3 Measurement Qubits. 
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1 
Let's take an example of the case when A = . where 


A is arational number with denominator as power of 
2. 


0. = (2e)2° — Yo2*/2°) 
= Cars = (o2" + Vil )/2") 
05 = Cag i = (792° - y12" + y227)/2°) 


Table 2 shows the probabilities, ||c(y)|| ? 
along with 09, 6;, 42 . 


Table 2 Probability of Measuring |y) in OPE Circuit for 
a Gate with Phase A = = Using 3 Measurement Qubits. 


[= [a [fom 


As the Table 2 illustrates, in the case when 


A= = for some integer k, 0; for c(k) is represented 
as in the following. 


t-1-i 
6, =| 2'k/2* - ». kao fe 
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t t-1-i 
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l # step 1: 


y 
= m +t t 
= Tl ; 2 Kplaeg] 2 
m'=0 
i 
m’' 
= 1 2 Kyglaped 
m!'=0 


It is shown that all 0; s of c(k) are integer 
multiples of m for some integer k when A= = ; 


Since cos* 6; = 1, the probability of measuring |k) 
is 1, whereas the probability of measuring |k’) for 
k'#k is 0. This result, shown in Table 2, is 
confirmed in the graph of Figure 2. 


4. EVALUATION 


The open-source quantum computing 
frameworks, such as IBM Qiskit [10] Google Cirq, 
Microsoft Q# and Rigetti Forest Quantum, has 
provided platforms for developing, simulating, and 
executing quantum algorithms [15]. 

We implement the QPE circuit using IBM 
Qiskit library as shown in Figure 4. The Qiskit 
circuit library provides OFT and PhaseEstimation 
classes implementing Quantum Fourier Transform 
and Quantum Phase Estimation. In this work, in 
order to illustrate the circuit composition, QPE 
circuit is composed of H, U gates and QFT class. 


The Python code in Figure 4 shows the source 
code used to build and simulate the QPE circuit 
using Qiskit. IBM released a major update of Qiskit 
from v0.46.0 to v1.0.1 as of February 2024. The 
following code is written in the new version v1.0.1. 


The circuit diagram of a quantum circuit built in 
the step 1 of Figure 4, is shown in Figure 5, 
generated using QuantumCircuit.draw() function. 


The simulation result in step 3 of code in Figure 
4 follows the probability distribution in Figure 3. 
The probability distribution confirms _ the 
computation in Table 1. 


Modifying the variable ¢ of code in Figure 4 
from 3 to 4, we built a QPE circuit for U with phase 


sf ; ee 
a for t=4 measurement qubits, which increases 


precision in the measurement. The simulation 
exhibits the probability distribution in Figure 6. 


# build a phase estimation circuit with 'phase' 

import numpy as np 

from qiskit import QuantumCircuit, \ 
QuantumRegister, ClassicalRegister 

from giskit.circuit.library import QFT 


t=3 

n=1 

qx = QuantumRegister(t,'qx’') 
qy = QuantumRegister(n, 'qy') 

c = ClassicalRegister(t,'c’) 

qpe = QuantumCircuit(qy, qx, c) 


phase = 1./3 
qpe.x(qy) 
qpe.barrier() 
qpe.h(qx) 
Un=1 
for qx_i in range(t): 

for i in range(Un): 

# controlled-Phase-Shift gate 
qpe.cp(2*np.pi* phase, qx[qx_i], qy); 
Un *=2 


# eigenstate |1> 


qpe.compose(QFT(t,inverse=True), \ 
range(n,n+t), inplace=True) \ 
qpe.measure(range(n,n+t),range(t))) 
# depiction of a quantum circuit 


qpe.draw(fold=-1, output='mpl') 


# step 2: 

# simulate the circuit as many as 'shots' times 

from qiskit.primitives.sampler import Sampler 

sampler = Sampler() 

job = sampler.run(qpe, shots=100000) 

result = job.result() 

prob = \ 
result.quasi_dists[0].binary_probabilities() 

print(prob) 


# step 3: 
# visualize the probability of measurement 
from qiskit.visualization \ 

import plot_distribution 
plot_distribution(prob) 


Figure 4 Python Source Code simulating OPE Circuit 
1 

for a Gate with Phase A = : Using 3 Measurement 

Qubits. 
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P (2n/3) 
qXo : 


P (2n/3) P (2n/3) 


Figure 5 Qiskit drawing of OPE circuit with n=1, t=3 


Table 3 Probability of Measuring y in QPE circuit for OPE Circuit for a Gate with Phase A = - Using 4 
Measurement Qubits. 


Computed Simulated 
pb |e fa [ae | aw Probability | Probability 
1 Z 4 8 
= = - = 0.004 0.004 
13 13 i3 13 
0001 — — — — 0.005 0.005 
8 


jto11) | on 0.004 0.004 
48 24” ad 
0.003 0.003 
0.003 0.003 
13 
|1111) = 0.003 0.003 
a" - Sn a eo 
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Quasi-probability 


S 

8 S 9S 
Figure 6 Qiskit Simulation Result of OPE Circuit for a 
Gate with Phase A = - USING 4 Measurement Qubits. 


The probabilities are shown in the last 
column in Table 3, labeled as simulated probability. 
Table 3 confirms that the result of the simulation is 
consistent with the computed probability from our 
analysis. 


1. 5 ane 
Because 3/8 close to ie the probability of 


measuring |y) = |0101) is highest and the 
probability decreases as the difference between y 
and 2° Ad increases, as it conforms to our intuition. 


4. EVALUATION 


Quantum phase estimation (QPE) is a vital 
component employed in many important quantum 
algorithms [16][17][18][19]. The ability to 
accurately estimate phase parameters is crucial for 
enabling quantum-enhanced computation in diverse 
fields ranging from cryptography [20][21][22][23] 
to material science [24][25][26]. 


In this paper, we have undertaken a comprehensive 
exploration of quantum phase estimation (QPE), 
focusing particularly on scenarios where the phase is 
not exactly represented as a rational number with a 
denominator in a power of 2. Through a rigorous 
mathematical analysis and precise computation of 
the quantum state of QPE in such cases, we have 
extended our understanding of this fundamental 
quantum algorithm. 


Our analysis contributes to a deeper appreciation of 
the behavior of QPE and its applications. By 
leveraging insights from number theory and 
trigonometry, we have derived explicit expressions 
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for the quantum state of QPE, providing valuable 
theoretical foundations for its practical 
implementation and optimization. In quantum 
chemistry and quantum simulation, where the phases 
represent energy eigenvalues or dynamics of 
quantum systems, exact or approximate computation 
of QPE enables accurate predictions of molecular 
properties, reaction rates, and material behavior. 
Similarly, in quantum machine learning and 
optimization, precise estimation of phases 
contributes to the development of efficient quantum 
algorithms for tasks like database search, pattern 
recognition, and optimization problems. 


The validation of our theoretical results through 
simulations using the Qiskit framework underscores 
the accuracy and efficacy of our formulations. The 
comparison between theoretical predictions and 
Qiskit simulations serves as a testament to the 
reliability and utility of our approach, offering 
practitioners and researchers a reliable mathematical 
framework for analyzing QPE. 


Moreover, our work contributes to the broader 
landscape of quantum computing _ research, 
highlighting the importance of understanding the 
behavior of fundamental quantum algorithms in 
various contexts. As quantum computing continues 
to evolve, with advancements in hardware, software, 
and algorithmic techniques, our insights into QPE 
pave the way for further innovation and exploration 
in the field. 


In conclusion, this paper advances our understanding 
of quantum phase estimation and similar quantum 
operations, offering mathematical insights and 
rigorous formulation of quantum states. 
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